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ABSTRACT 

Context. We present the results from Suzaku observations of the merging cluster of galaxies CIZA J2242.8+5301 at z=0.192. 

Aims. To study the physics of gas heating and particle acceleration in cluster mergers, we investigated the X-ray emission from CIZA 
J2242.8+5301, which hosts two giant radio relics in the northem/southern part of the cluster. 

Methods. We analyzed data from three-pointed Suzaku observations of CIZA J2242.8+5301 to derive the temperature distribution in 
four different directions. 

Results. The Intra-Cluster Medium (ICM) temperature shows a remarkable drop from 8.5 4 jj' 8 keV to 2.7 4 ]J J keV across the northern 
radio relic. The temperature drop is consistent with a Mach number Ai„ = 2.7 4 !H and a shock velocity v^ock-.n = 2300 4 ^jj km sA 
We also confirm the temperature drop across the southern radio relic. However, the ICM temperature beyond this relic is much higher 
than beyond the northern one, which gives a Mach number M s = l-7 4 jj' 4 and shock velocity v s hock.s = 2040km sA These results 
agree with other systems showing a relationship between the radio relics and shock fronts which are induced by merging activity. We 
compare the X-ray derived Mach numbers with the radio derived Mach numbers from the radio spectral index under the assumption 
of diffusive shock acceleration in the linear test particle regime. For the northern radio relic, the Mach numbers derived from X- 
ray and radio observations agree with each other. Based on the shock velocities, we estimate that CIZA J2242.8+5301 is observed 
approximately 0.6 Gyr after core passage. The magnetic field pressure at the northern relic is estimated to be 9% of the thermal 
pressure. 

Keywords, galaxies: clusters: individual (CIZA J2242.8+5301)—galaxies: intergalactic medium — shockwaves — X-rays: galaxies: 
clusters 


1. Introduction 


Based on the framework of the hierarchical structure formation 
theory, clusters of galaxies grow through merging events with 
smaller objects and accretion flows from large scale structure 
filaments (Voit 2005; Kravtsov & Borgani 2012). In particu¬ 
lar, merging events between massive clusters of galaxies are the 
most energetic events in the Universe since the Big Bang, releas¬ 
ing as much as ~10 64 erg (Markevitch & Vikhlinin 2007). This 
gravitational energy is converted into heating of the Intra-Cluster 
Medium (ICM) and particle acceleration by shocks likely via the 
mechanism of diffusive shock acceleration (DSA: Blandford & 
Eichler 1987). 


Radio observations show the existence of diffuse non- 
thermal emission in clusters such as radio halos, relics and mini¬ 
halos (for a review, see Feretti et al. 2012; Brunetti & Jones 
2014, and references therein). Among these radio structures, 
radio relics show remarkable features such as a diffuse (Mpc- 
scale), elongated shape, and they are typically located at cluster 


outskirts. In addition, radio relics are usually seen in merging 
clusters. Therefore, radio relics are considered to be the result 
of synchrotron emission from relativistic electrons which are ac¬ 
celerated by shocks induced by cluster mergers (Ensslin et al. 
1998). 

Because shocks play a fundamental role in cluster evolution, 
it is crucial to understand their properties (Mach number, shock 
velocity, shock acceleration efficiency, etc.). In particular, at 
high Mach numbers (Al ~ 10 3 - 10 4 ) DSA is known to have 
enough efficiency to accelerate particles from a thermal distribu¬ 
tion to the high energy regime. This is confirmed sby observa¬ 
tions of several supernova remnants (e.g., Koyama et al. 1995). 
However, shock fronts induced by cluster mergers are expected 
to have a much lower Mach number (AT ~ 2-5). The accel¬ 
eration efficiency of these low-Mach number shocks is thought 
to be too low to reproduce the observed radio brightness (Kang 
et al. 2012). However, recent work by Guo et al. (2014a,b), us¬ 
ing detailed particle in cell simulations indicates, that electron 
acceleration can be efficient at low-Mach number shocks in the 
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Fig. 1 . Left: Suzaku observations of CIZA J2242.8+5301 overlaid on a 1.4 GHz radio intensity map in color scale (van Weeren et al. 2010). 
The Suzaku XIS fields of view are shown as white squares. Right: X-ray image of CIZA J2242.8+5301 in the energy band 0.5-2.0 keV, after 
subtraction of the NXB with no vignetting correction and after smoothing by a 2-dimensional gaussian with <x =8 pixel =8.5". The 1.4 GHz radio 
emission is shown with yellow contours. 


ICM. More detailed studies of cluster shocks are highly desired 
for understanding their roles during merging activity. 

Most evidence for shock fronts in clusters of galaxies is 
based on X-ray imaging studies (for a review, see Markevitch 
2010, and references therein). Up to now, clear evidence for 
shock fronts are reported in IE 0657-56 (Bullet cluster: Marke¬ 
vitch et al. 2002), A520 (Markevitch et al. 2005), A754 (Macario 
et al. 2011), A2146 (Russell et al. 2010, 2012), A521 (Bourdin 
et al. 2013) and A2034 (Owers et al. 2014). 

X-ray follow-up observations of the northwest radio relic 
in Abell 3667 revealed that the ICM temperature and surface 
brightness distribution drops across the relic, which indicates 
the existence of a shock front across the relic (Finoguenov et al. 
2010). The correspondence between radio relics and shocks is 
also confirmed in other clusters (Macario et al. 2011; Akamatsu 
et al. 2012a,b; Ogrean & Briiggen 2013; Bourdin et al. 2013; 
Akamatsu et al. 2013). Using the Rankine-Hugoniot jump con¬ 
ditions, the estimated Mach number from the temperature jump 
typically spans M ~ 1.5 - 3.0 (Akamatsu & Kawahara 2013). 
These values are consistent with prediction from hydrodynami- 
cal simulations (Miniati et al. 2000; Ryu et al. 2003). 

This paper reports the results of deep X-ray observations of 
the cluster CIZA J2242.8+5301 with Suzaku. It has been dis¬ 
covered in the Clusters in the Zone of Avoidance (CIZA) survey 
(z = 0.192: Kocevski et al. 2007). CIZA J2242.8+5301 has a 
clear double radio relic at a scale of several Mpc and with a 
remarkably narrow width of 50 kpc for the northern relic (van 
Weeren et al. 2010). The relic is strongly polarized at the 50 
to 60% level, indicating a well-ordered magnetic field, which is 
well aligned with the shock plane. 

The northern relic shows an injection index at the edge of 
the relic of -0.6+0.05 (van Weeren et al. 2010; Stroe et al. 2013), 
which corresponds to a Mach number of M ~ 4.6. The spectral 
index for the south relic is also reported as -1.29 + 0.05, which 


corresponds to a Mach number of Al ~ 2.8 (Stroe et al. 2013). 
Stroe et al. (2014b) also report the detection of radio emission 
at the high frequency band (~ 16 GHz) and pose some questions 
for the diffusive shock acceleration predictions of the ageing ef¬ 
fect on relativistic electrons. A recent deep Chandra observation 
of CIZA J2242.8+5301 revealed several edges in the X-ray sur¬ 
face brightness distribution. Interestingly, the edges are found in 
the downstream regions of the shocks assumed to be associated 
with the radio relics, and have no counterpart in the temperature 
distribution (Ogrean et al. 2014). 

Previous Suzaku observations of CIZA J2242.8+5301 
showed a remarkable temperature jump from 8.3 keV to 2.1 keV 
across the northern relic (Akamatsu & Kawahara 2013). The 
Mach number estimated from the Rankine-Hugoniot condition 
is Mx = 3.2 + 0.5. On the other hand, as described above, van 
Weeren et al. (2010) reported a somewhat higher Mach number 
A^Radio = 4.6^^. Similar higher values, Af ra dio > Mx, have 
been found for two other relics (1RXS J0603: van Weeren et al. 
2012; Ogrean et al. 2013, Itahana in prep). Recently, Sarazin 
et al. (2014) reported results of deep XMM-Newton observations 
of the northeast radio relic in Abell 3667. They derive a Mach 
number of 2.09+0.09, which is clearly much smaller than the 
Mach number expected from the radio spectral index (Spectrum 
index q'=- 0.6+0.2, M=4.6+0.3: Hindson et al. 2014). 

However, more recently Stroe et al. (2014a) reported the 
possibility of underestimating the radio spectral index from 
spatially-resolved spectral fitting analysis for the northern radio 
relic in CIZA J2242.8+5301 and determine a new spectral index 
of 0.77 ^qq 2- The Mach number estimated from the new radio 
spectral index is M - 2.9 '[[[l!, which is now in agreement with 
X-ray observations (Akamatsu & Kawahara 2013; Ogrean et al. 
2014). 

The combination of X-ray and radio observations gives us 
the opportunity to solve many problems related to shock physics, 


Article number, page 2 of 12 










H. Akamatsu et al.: S uzaku X-ray observations of CIZA J2242.8+5301 


Table 1. Suzaku observations log of the CIZA2242 


Region 

(R.A, DEC) 

Exposure 

(ks) 

North 

(340.74,53.16) 

82.2 

East 

(340.75, 52.92) 

74.0 

West 

(340.51,52.95) 

77.4 

OFFSET 

(339.29,52.67) 

38.2 


in particular the nature of the shock fronts themselves and their 
role in the acceleration and heating processes in the ICM. From 
this point of view, CIZA J2242.8+5301 is an important target to 
study at X-ray wavelengths. 

We use Hq — 70 km s -1 Mpc -1 , Qm = 0.27 and Qa = 0.73, 
which gives a scale of 190 kpc per arcminute at z - 0.192. We 
employ solar abundances defined by Lodders (2003) and Galac¬ 
tic absorption with N H: totai = 46.0 x 10 20 cm~ 1 2 (Willingale et al. 
2013). Unless otherwise stated, the errors correspond to 68% 
confidence for a single parameter. 


2. Suzaku observations & Spectrum analysis 

Suzaku carried out three observations of CIZA J2242.8+5301 
(Figure 1 left) and one observation of an OFFSET region, dur¬ 
ing the Suzaku AO-6 (PI: Kawahara) and AO-8 (PI:van Weeren) 
phase. We refer to these pointings as North, East and South, re¬ 
spectively. They cover the whole cluster beyond the radio relics. 
The observation log is given in Table 1. All observations were 
performed with the normal 5x5 and 3x3 clocking mode 1 (with¬ 
out burst or windows options). 

The XIS instrument consists of 4 CCD chips: one back- 
illuminated (BI: XIS1) and three front-illuminated (FI: XIS0, 
XIS2, XIS3) chips (Mitsuda et al. 2007). XIS2 was not oper¬ 
ational due to damage from a meteoroid strike 2 . Data reduc¬ 
tion was done with HEASOFT version 6.15 and CALDB version 
20140203. 

The XIS event lists created by the rev 2.5 pipeline processing 
were filtered using the following additional criteria: a geomag¬ 
netic cosmic-ray cut-off rigidity (COR) > 8 GV, and an eleva¬ 
tion angle ELV > 12°. We applied additional processing for the 
XIS 1 detector to reduce the non-X-ray background (NXB) level, 
which has increased after the change in the amount of charge 
injection from 2 keV to 6 keV in lune 2011. That operation 
improved the high-energy response with a negligible loss in the 
low-energy performance 3 . The detailed processing procedures 
are the same as those described in XIS analysis topics 4 . The 
5x5 and 3x3 editing modes data formats were added. 

The observed spectrum was assumed to consist of optically- 
thin thermal plasma emission from the ICM, and emission due 
to the Galactic foreground. Cosmic X-ray background (CXB) 
and NXB. In order to investigate the properties of the clus¬ 
ter plasma outside of the radio relics, an accurate estimation of 
the foreground-background emission is essential. Here we use 
a Suzaku offset observation, which is located 1 degree east of 
CIZA 12242.8+5301, to estimate the sky background spectra. 


1 http://www.astro.isas.ac.jp/suzaku/doc/suzaku_td/node 10.html 

2 http://www.astro.isas.ac.jp/suzaku/doc/suzakumemo/suzakumemo- 
2007-08.pdf 

3 http://www.astro.isas.ac.jp/suzaku/doc/suzakumemo/suzakumemo- 
2010-07v4.pdf 

4 http://www.astro.isas.jaxa.jp/suzaku/analysis/xis/xisl_ci_6_nxb/ 


Sky background spectrum in CIZA2242 OFFSET region 
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Fig. 2. The spectrum of the OFFSET field used for the background 
estimation, after subtraction of the NXB and the point sources. The XIS 
BI (Black) and FI (Red) spectra are fitted with CXB + Galactic com¬ 
ponents (LHB. MWH and HF) ( apec+ phabs(apec+powerlaw+apec)). 
The CXB spectrum is shown with a black curve, and the LHB, MWH 
and HF components are indicated by green, blue and gray curves, re¬ 
spectively. 

In the fitting procedure for the cluster emission, we com¬ 
bined the spectra from the FI detector (XISO, 3) and BI and FI 
spectra were fitted simultaneously. The NXB subtracted spec¬ 
trum is fitted with a model of the ICM and the sky-background 
components. Details about background estimation and the ICM 
components across the northern/southern radio relics will be de¬ 
scribed in subsections 2.1, 2.3 and 2.4. 

2.1. Estimation of the background spectra 

The sky X-ray background is divided into four components: un¬ 
absorbed thermal emission from the Local Hot Bubble (LHB: 
kT ~0.1 keV), absorbed thermal emission form the Milky Way 
Halo (MWH: kT ~ 0.3 keV), the Hot Foreground (HF: kT ~ 0.6 
keV) and the CXB. 

In order to determine the level of the sky X-ray back¬ 
ground, we first analyzed a Suzaku offset observation of CIZA 
J2242.8+5301 (obsID:806002010). We assume that these sky 
background components are distributed uniformly over the 
Suzaku FOV. We used uniform ancillary response files (ARF) 
for the sky background components. The ARF for the sky back¬ 
ground was generated by using xissimarfgen (Ishisaki et al. 
2007) by assuming uniform-sky emission over a circular region 
with a radius of 20'. Hereafter, we call this uniform-ARF. 

The NXB component was reconstructed from the XIS Night 
Earth database with xisnxbgen in FT00LS (Tawa et al. 2008). We 
accumulated the data for the same detector area and the same dis¬ 
tribution of COR2 as the observations. To increase the signal to 
noise ratio, we applied data selection with COR2 > 8 GV (see 
Figure 3 in Tawa et al. 2008). To adjust for the long-term vari¬ 
ation of the XIS background due to radiation damage, we select 
the night Earth data within 300 days before and after the period 
of the observation. Using the wavdetect function in the CIAO 
package 5 , we searched for point sources in the Suzaku image. 
Point sources were detected down to ~ 1x10 13 erg cm -2 s 1 
in the 2-10 keV band. We exclude the point sources within a 
1-2 arcmin radius in order to take account for the point spread 
function (PSF) of the Suzaku XRT (Serlemitsos et al. 2007). 

5 http://cxc.harvard.edu/ciao/ahelp/wavdetect.html 
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Table 2. Best-fit background parameters 



Local Hot Bubble 

Milky Way Halo 

Hot Foreground 

CXB 



kT (keV) 

norm d 

[xlO -3 ] 

kT (keV) 

norm a 
[XlO- 3 ] 

kT (keV) 

norm a 
[XlO” 3 ] 

r 

norm? 

[xlO- 4 ] 

X 2 /d.o.f 

Casel 

0.08 (fix) 

24.2 ± 2.2 

O 

o 

1 + 

o c 

b c 

— k; 

4.9 + 0.4 

- 

- 

1.41 (fix) 

9.4 + 0.4 

166/133 

Case 2 

0.08 (fix) 

43.2 + 5.0 

0 27+ +0 01 
-0.04 

6.1 + 1.5 

'nI- oc 

o c 

d d 

+ l 

OO 

VO 

d 

1 o +0 - 6 

_ -0.3 

1.41 (fix) 

9.2+ 0.4 

159/132 


a: Normalization of the apec component scaled by a factor 20 x 20;r = 400;r (see text). 

norm=-^~ f n e n H dV/(4n(l +z 2 )D\)x 10 14 cm 5 arcmin 2 , where D A is the angular diameter distance to the source. 

b: The CXB intensity normalization in Kushino et al. 2002 is 9.6 x 10 -4 for T = 1.41 in units of photons keW 1 cirr 2 s _1 at 1 keV. 


o 

o 

o 

CM 


O 

LU 

0 
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CIZA J2242.8+5301 merging direction 



Fig. 3. Left: X-ray image of the CIZA J2242.8+5301 in the energy band 0.5-5.0 keV, after subtraction of the NXB with no vignetting correction 
and after smoothing by a 2-dimensional gaussian with cr = 8 pixel =8.5". The identified point sources by XMM-Newton are highlighted with 
the green circles (see text). The 1.4 GHz radio emission is shown with yellow contours. The dashed cyan circle indicate the center which was 
determined by visually fitting a circle to the radio relics. Right: Temperature profile of CIZA J2242.8+5301 from the south to north direction. The 
dotted gray lines show WSRT 1.4 GHz radio emission. 


For the spectral fits, we used the XSPEC version 12.8.0 pack¬ 
age 6 . The model for the sky-background components is de¬ 
scribed as ApecLUB+phabs(ApecM'WH+ApeciiF+powerlawcxB ) 
in XSPEC representation. In this model, Apec and phabs repre¬ 
sents thin thermal plasma emission model (Smith et al. 2001) 
and Galactic absorption toward the target, respectively. We fixed 
the abundance and the redshift of each thermal component to 1.0 
and 0.0, respectively. Further we fixed the temperature of the 
LHB and photon index of the CXB components to 0.08 keV and 
r=1.41 (Kushino et al. 2002), respectively. We used the spectra 
in the 0.5-7 keV range for the BI detectors and 0.8-7 keV for 
the FI detector. The result is consistent with the CXB intensity 
with the Kushino et al. (2002) level and the MWH tempera¬ 
ture of 0.27 /q^ keV is consistent with typical values in other 
fields (Yoshino et al. 2009). A fit with the HF component fixed 
to zero (easel) is sightly poorer than a fit with its temperature 
and normalisation as free parameters (case2). The resulting fit 
parameters are shown in Table 2. Basically these models repro¬ 
duce all features of the observed spectrum well, CXB intensity, 
low energy (< 1 keV) sky-background. We have checked how 
the HF component improves the fitting with an F-test, and ob- 

6 http://heasarc.nasa.gov/xanadu/xspec/ 


tained the probability value of the F-test ~ 0.001. Thus, we em¬ 
ployed the parameters of case 2 in Table 2. The best-fit spectrum 
model (case2) is shown in Figure 2. 

2.2. ICM emission along the merging direction 

In order to investigate the cluster temperature distribution in the 
merging direction (from north to south), we divided the cluster 
into 11 box regions as shown in Figure 3. For the ICM emission 
around the radio relic, we employed an absorbed single temper¬ 
ature thermal plasma model (phabs x apec). In all regions, we 
fixed the metal abundance and the redshift to 0.3 (Fujita et al. 
2008; Werner et al. 2013) and 0.192 (Kocevski et al. 2007), re¬ 
spectively. Because the ICM emission has a different spatial dis¬ 
tribution from the sky-background, we generated an ARF by us¬ 
ing the Suzaku XIS image (0.5-8.0 keV band) as input surface 
distribution for xissimarfgen. Based on Sec. 2.1, we modeled 
the sky-background emission using the uniform ARF. Employ¬ 
ing multiple ARFs is enabled with XSPEC version 12. The nor¬ 
malization of the LHB component, the normalization and tem¬ 
perature of the MWH component, and the normalization of the 
power-law model for the CXB component were allowed to vary 
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Outside the northern radio relic 



Energy (keV) 


Inside the northern radio relic 



Energy (keV) 


Fig. 4. NXB subtracted spectra in each annular region. The XIS BI (Black) and FI (Red) spectra are fitted with the ICM model (phabs x apec), 
along with the sum of the CXB and the Galactic emission of Case 2 (apec + phabs(apec + apec + powerlaw)). The ICM component is shown 
with a magenta line. The CXB component is shown with a black curve, and the LHB, MWH and HF emissions are indicated by green, blue and 
gray curves, respectively. The sum of all sky background components is shown by the orange curve. The Left panel shows the spectrum of the 
outside region. Right panel shows the inside region. 



24 22:44:00 43:36 43:12 48 24 42:00 

RA (J2000) 



0.23 0.26 0.34 0.48 0.67 0.92 1.2 1.6 2 2.5 3 



Fig. 5. Left: Suzaku XIS image for the northern part of CIZA J2242.8+5301. Right: Radial profile of the ICM temperature for the north 
direction. The present best-fit values with statistical errors are shown with black crosses. Gray dotted crosses indicate previous Suzaku results, 
which do not include HF component (Akamatsu & Kawahara 2013). The range of uncertainties due to the combined 3% variation of the NXB 
level and the maximum/minimum fluctuation in the CXB is shown by two green dashed lines.The black dotted lines show the WSRT 1.4 GHz 
radio intensity scaled arbitrary in flux. The higher peak corresponds the northern radio relic and the second peak is a point source just behind the 
radio relic. 


within the range of the errors in the background estimation (ta¬ 
ble 2). We carried out spectral fits to the pulse-height spectra in 
each annular region separately. In the simultaneous fit of the BI 
and FI data, only the normalizations were allowed to be differ¬ 
ent between them, although we found that the derived normal¬ 
izations were consistent within 15% The resulting temperature 


profile is shown in Figure 3. The temperature profile shows a 
significant drop at r --6 and 6 arcmin, respectively. The tem¬ 
perature drop strongly indicates the presence of shock fronts at 
those locations. Here we note that the temperature profile was 
smeared by the point spread function of the Suzaku XRT (half 
power diameter ~ 2 arcmin: Serlemitsos et al. 2007). Therefore, 
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we can not resolve structures of scales less than an arcmin. In the 
next section, we will investigate in detail the properties around 
the radio relics. 

2.3. ICM emission across the northern radio relic 

To investigate the detailed ICM properties associated with the 
radio relics we extracted pulse-height spectra in two annular re¬ 
gions whose boundary radii were 4.0'-7.0' and 7.5'-11.5' for 
the region inside and outside the radio relic with the center at 
(22:42:41.9, 53:03:00.0), shown in Figure 5. To reduce contam¬ 
ination from the brighter region, we take the region outside the 
radio relic 0.5' away from the outer edge of the inner region. 
To maximize the signal to noise ratio outside the northern relic 
region, we used the spectra in the 0.5-4 keV range for the BI 
detector and 0.8-4 keV for the FI detectors. 

To investigate a possible systematic error in the tempera¬ 
ture estimation, we consider a 3% fluctuation of the NXB level 
and fluctuations of the CXB intensity as shown below. It is 
important to exclude the contribution of point sources in the 
region of interest, to reduce the systematic error of the CXB 
caused by fluctuations of unresolved sources. For the elimina¬ 
tion of the point sources, we used the archival data of XMM- 
Newton (OBSID=0654030101). As shown by small (green) 
circles in Figure 3, we detected several point sources and sub¬ 
tracted 10 point sources with 2.0-10.0 keV fluxes higher than 
S c = 1 x 10 14 erg cm -1 s 1 . According to previous studies 
of the CXB with the ASCA satellite (Kushino et al. 2002), the 
CXB flux fluctuation in the field of view of the ASCA GIS (0.5 
deg 2 ) is 6.5% in the 2-10 keV band. We scaled the measured 
fluctuations from our flux limit and the observed area follow¬ 
ing Hoshino et al. (2010). The estimated fluctuations are 32% 
and 18% for inside and outside of the radio relic regions, re¬ 
spectively. We repeated all the spectral fits by fixing the CXB 
intensity at the upper and lower boundary values. 

The pulse-height spectra and the best-fit models for both an¬ 
nular regions are shown in Figure 4. We obtained fairly good fits 
for all the regions with reduced^- 2 values less than 1.2. The re¬ 
sultant ICM temperature with different background models are 
listed in table 5. The temperature changes significantly from 
fcTinside = 8.5^°;g keV to AToutside = 2.7 + _%' A keV. Even though we 
use different sky background models, these results are consis¬ 
tent with previous Suzciku results (/f7’ mS ide = 8.3 ± 0.8 keV and 
Ar ou t s ide = 2.1 ± 0.4 keV: Akamatsu & Kawahara 2013). These 
ICM temperatures are consistent with the ones estimated from 
the box regions shown in Figure 3. 

For the inside and outside the northern relic, box region 
shows A7i ns ide = 9.7 +1.3 keV and Ar outs id e = 3.4 ±1.1 keV, re¬ 
spectively. The ICM temperature of the post-shock region also 
agrees with a previous XMM-Newton observation (A7) nsK i c ~ 
8-11 keV: Ogrean et al. 2013). Figure 5 shows the radial profile 
of the ICM temperature across the northern radio relic. Black 
crosses show the best fit ICM temperature of each annular re¬ 
gion. Grey dotted crosses represent previous results (Akamatsu 
& Kawahara 2013), which do not include the HF component. 

One can also estimate the shock compression from the jump 
in the X-ray surface brightness across the relic. The surface 
brightness was derived from the observed spectra. We checked 
the surface brightness across the well defined northern radio 
relic. The resulting un-absorbed X-ray surface brightness (0.5- 
2.0 keV) shows a factor 8 difference across the relic, which cor¬ 
responds to a shock compression of C ~ V8.22 = 2.86, includ¬ 
ing the difference in emissivity for the two regions. This value 


Table 3. Best-fit ICM temperature across the northern radio relic with 
different background model (unit of keV) 



Casel 

Case2“ 

Case3 b 

Outside N relic 

9 7 + 0.5 

-Q.4 

9 7+0./+0.5 

^ '-Q.4-Q.3 
o c+0.8+0.2 
-0-6—0.2 

2 6 +uy 

-0.5 

Inside N relic 

0 ° 

u> 

1 1 

3 C 

8 .6+0.9 


a: The first error is the statistical one and the second one 
is due to systematics one. 

b: Case 2 background model with > 1 keV energy band. 


Table 4. Best-fit ICM temperature across the southern radio relic with 
different background model (unit of keV) 



Casel 

Case2 a 

Case3 b 

Beyond S relic 

4 4+ () -8 

-Q-3 

c rt+O.y + l.l 

-0.6-0.9 

5 2 +2 - u 

J -0.9 

Inside S relic 

o q+0.4 
° -0.3 

9.0+0.5+0.4 

9.0+0 .6 


a: Idem as above, 
b: Idem as above. 


matches the compression parameter (C ~ V7.41 = 2.72) derived 
across the rectangular sector, which is shown Figure 3. These 
surface brightness jumps are consistent with the shock compres¬ 
sion predicted based on the temperature jump (Table.5: C ~ 
2.83^q 2 q). This agrees with the results for shocks found on other 
clusters (for instance. Bullet: Markevitch et al. 2002, Shimwell 
et al. 2015, A520: Markevitch et al. 2005, A754:Macario et 
al. 2011, RXJ1314.4-2515 Mazzotta et al. 2011). However, we 
note that our estimation based on our X-ray spectroscopic anal¬ 
ysis represent the upper limit to the density jumps because the 
surface brightness profiles of clusters of galaxies typically have 
steep power-law profiles. It is difficult to make concrete state¬ 
ments on this issue due to the limited spatial resolution of the 
Suzaku XRT. 

2.4. ICM emission around the southern radio relic 

Next we investigate the ICM properties around the southern ra¬ 
dio relic with our new Suzaku observations. To increase the sig¬ 
nal with our limited observation time, the West and East obser¬ 
vations overlap at the southern radio relic. We analyzed these 
data in the same way as the northern one. Similar as for the 
northern relic, we take a one arcmin gap to avoid contamina¬ 
tion from the brighter part of the cluster due to the PSF of the 
Suzaku XRT. To investigate possible structures in the W and E 
directions, we also evaluate the ICM temperatures by shifting 
the annuli by 0.5'. The resulting ICM temperatures are shown 
by the gray crosses in the Figure 7. 

We used the spectra in the 0.5-7 keV range for the BI detec¬ 
tors and 0.8-7 keV for the FI detector. For point source exclu¬ 
sion, we only used Suzaku data because XMM-Newton/Chandra 
observations did not completely cover the full Suzaku FOV. We 
excluded point sources by using CIAO to a flux detection limit 
of S c — lx 10~ 13 erg cmT 2 s -1 . Based on the flux detection limit, 
we estimate the possible fluctuation of the CXB intensity in each 
region. The estimated fluctuations span 35-60 %. 

We successfully obtained the ICM temperature across and 
beyond the south relic for the first time. The resulting ICM tem¬ 
perature across the radio relic region is shown in Table 4. The 
best-fit temperature profiles are shown in Figure 7. Toward the 
southern radio relic, the ICM temperature show slightly increas¬ 
ing profile from 8.0+0.3 keV to 9.0+0.5 keV and shows a signif¬ 
icant drop to 5.0^gg keV. Similar to the results for the northern 
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Fig. 6 . Radial profile of the ICM temperature for the southern direction. Colors are the same as in Figure. 5. 
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Fig. 7. Radial profile of the ICM temperature for the South (left), East and West direction. The Black crosses show the best-fitted values and the 
gray crosses show the the results of slightly offset regions. The gray dotted curves in the right panels show the temperature profile (Bums et al. 
2010) (see text). In the left panel, we removed the radio point sousrce located just behind the relic. 


shocks/relic, these temperatures are consistent with those de¬ 
rived for in the rectangular sector (Figure 3), although the sta¬ 
tistical errors are large. This might be due to the fact that this 
region is located too far from the radio relic, the bin is placed 
about 760 kpc from the peak of radio emission. Given the ex¬ 
pected decrease in temperature with radius this would explain 
the lower temperature measured there. Combined with the fluc¬ 
tuation of the NXB level (3 %: Tawa et al. 2008), we investigated 
the systematic error. The results are shown by green dotted lines 
in Figure 7. Due to the high point source flux detection limit, the 
fluctuations of the CXB intensity are large. The emission inten¬ 
sity outside of the southern radio relic is almost comparable to 
the sky background (Figure 6). Therefore, the systematic errors 
outside of the southern radio relic are larger than the statistical 
error. 


3. Discussion 

We performed deep Suzaku observations of the double radio relic 
cluster CIZA J2242.8+5301. We measure the ICM temperature 
profiles beyond the relics were measured. The profiles show a 
significant drop across the radio relics, which implies the pres¬ 
ence of shock fronts. We discuss the full cluster temperature 
structure by comparing to other clusters. We evaluate the shock 
properties (Mach number, compression factor and propagation 
speed), the possibility of non-equilibrium effects due to shock 
heating and we compare these results with the radio observa¬ 
tions. 

3.1. Large scale ICM temperature structure 

As suggested by numerical simulations, clusters grow via merg¬ 
ing with sub-clusters. These merging activities have a huge im- 
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Table 5. Shock Properties at the northern and southern radio relic 



Ta 

1 2 

r r a 

1 1 

Mach No. b 

^ shock 

Compression 

Power-law slope 9 

Spectrum index 


(keV) 

(keV) 

M 

(km s 1 ) 

C 

P 

a 

Northern 

0 c+0.8 

-0.6 

9 '7+0.9 

- 0.5 

9 9+0.7 

^-0.4 

2300 +7O(J 

2 84+U+3 

-0.24 

2.64 + 0.30 

0.81+0.15 

Southern 

9.0 + 0.6 

5 1 + ^ 

1 7+0.4 
1,/ -0.3 

2040 +55 ° 

ZWU -41D 

2 Q 1+0.45 
-0.34 

3 96 + 1-32 
- 0.99 

1 48+0-66 
i -^°-0.49 


1,2 indicate pre-shock and post-shock, respectively. 

a: The error ranges calculated by considering statistical and systematic errors (see text in Sect.3.2 for details). 


b: From the Rankine-Hugoniot temperature jump condition ^ = iAdFdiW - 

c: The shock speed v s h oc k — At ■ v s , d: Expected Power-law slope p — (C + 2)/(C - 1) 



Fig. 8. Scaled projected temperature profiles compared with relaxed 
clusters (Gray: see Fig. 8 in Reiprich et al. 2013). The profiles were 
normalized to the average temperature. The r 2 oo value was derived 
from Henry et al. (2009). Blue, red, yellow and magenta crosses show 
the temperature profile of the North, South, East and West direction, 
respectively. 


pact on the temperature structure (see for a review Markevitch 
& Vikhlinin 2007). Recent Suzaku observations reveal the ICM 
temperature profiles up to the virial radius (see for a review 
Reiprich et al. 2013). However, these studies mainly focus on 
relaxed clusters and studies of merging cluster outskirts are still 
limited except for a few examples (e.g., Ibaraki et al. 2014). 
Therefore, the impact of merging activity on the cluster struc¬ 
ture remains poorly known. 

The new Suzaku observations of CIZA J2242.8+5301, in 
combination with Suzaku observations taken in 2011 (Akamatsu 
& Kawahara 2013), cover the full cluster region. As shown in 
Figure 3, CIZA J2242.8+5301 has a peculiar temperature pro¬ 
file along the merging axis (from north to south). Contrary to 
relaxed clusters (Reiprich et al. 2013), CIZA J2242.8+5301 has 
an almost flat radial profile with kT ~ 8 keV (Ogrean et al. 2013, 
and this work). In addition, across the radio relics, the tempera¬ 
ture profile also shows a significant drop. These observed tem¬ 
perature profiles strongly indicate the presence of a shock front 
and shock heating on the ICM. We discuss the properties of the 
shock fronts in the next subsection. 

We plotted a model temperature profile proposed by Burns 
et al. (2010) as gray dotted curves in the right panel of Figure. 7. 
Their model profile reproduces well the measured ICM tempera¬ 
ture structure of relaxed clusters of galaxies by the Suzaku satel¬ 
lite (Reiprich et al. 2013). For the calculation of the model pro¬ 
file, we adopt 7' avg = 8.0 keV and r 2 oo = 12 arcmin, respec¬ 


tively. Here, r 2 oo 7 is estimated from Henry et al. (2009). We also 
compare the measured temperature profiles with previously mea¬ 
sured clusters (Figure 8). The temperature profile in the West di¬ 
rection agrees with the model profile. On the other hand, the pro¬ 
file in the East direction shows an excess. The temperature ex¬ 
cess in the East direction is consistent with the previous XMM- 
Newton results (Ogrean et al. 2013). They report the presence 
of a high temperature (kT ~ 13 keV) region at the eastern part of 
the cluster. In our temperature profile, we do not find such high 
a temperature region. However, due to the limited PSF of the 
Suzaku XRT, we may miss such a component. In addition, such 
a high temperature (kT >10 keV) is beyond the effective energy 
band of XIS and of XMM-Newton. Therefore tight constraints 
on high temperature components are difficult to obtain for cur¬ 
rent X-ray CCD detectors. NuSTAR 8 9 (Harrison et al. 2010) and 
the upcoming ASTRO-H 9 (Takahashi et al. 2012) satellite will 
help to detect such high temperature components with their hard 
X-ray imaging and spectral capabilities. 

3.2. Shock properties 

The Suzaku temperature profiles of CIZA J2242.8+5301 show 
clear drops across the relics, indicating the presence of shock 
fronts. We evaluate the shock properties at both radio relics 
based on the Suzaku results. Hereafter we used the error in¬ 
cluding systematic errors due to the variation of the NXB level 
and the maximum/minimum fluctuation in the CXB defined as 

cr err = ^ Jo~s lat + 0"s y s- We cann °t determine a one-to-one rela¬ 
tionship between the edge of the radio relics and the shock fronts 
because of the limited PSF of the Suzaku XRT. 

The Mach number can be obtained by applying the Rankine- 
Hugoniot jump condition (Fandau & Fifshit’s 1959) 

T 2 5A+ + 14A+ - 3 

74 “ 16M 2 ’ ( } 

assuming the ratio of specific heats as y = 5/3. Here, 1,2 indi¬ 
cate pre-shock and post-shock, respectively. Table 5 shows the 
resultant Mach numbers based on the temperature jump, whose 
values are Af n = 2.7/°'^, A4 S = 1.7/;° 3, respectively (Table 5). 

Using the measured temperature in the pre shock region 
(^7’ pre:n ~ 2.7 keV, kT pie:s ~ 5.1 keV), the sound speed is 
c s:n ~ 850 km s -1 , c s:s ~ 1170 km s -1 . With these Mach num¬ 
bers, the expected shock compression C and the shock propaga¬ 
tion speed (v s hock = A1 • c s ) are estimated as C n = 2.84/°^, 

7 r 2 oo is the radius where the mean total density of the cluster is 200 
times the mass density of the Universe with critical density. 

8 http://www.nustar.caltech.edu/ 

9 http://astro-h.isas.jaxa.jp/index.html.en 
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Fig. 9. Smoothed 0.5-2.0 keV band X-ray image (red) and WSRT 1.4 
GHz image of CIZA J2242.8+5301 (cyan). The thin yellow lines depict 
the approximate locations of the shock fronts confirmed by Suzaku. 

C s = 2.01+°;£j and v shock:n = 2300+™° km s -1 , v shock:s = 
2040+4™ km s _1 , respectively. For the northern radio relic, the 
estimated shock speed is in agreement with the radio observa¬ 
tions ( v s hock-.n ~ 2500 km s _I : Stroe et al. 2014a). Here we esti¬ 
mate the shock compression C from the Mach numbers derived 
from the temperature drops. Although the Mach number is com¬ 
parable to that of the Bullet cluster (Al ~ 3.0), the shock velocity 
in CIZA J2242.8+5301 is much smaller than for the Bullet clus¬ 
ter (4500 km s -1 ) and is comparable to A520 (2300 km s _1 ) and 
other systems (e.g. A2034: 2057 km s _1 ). The relative gas flow 
speed to the shock front is v gaS :n = v s i l0 ck/C ~ 810 km s _1 and 
Vgas-.s ~ 1000 km s _1 . The estimated downstream velocity for the 
northern shock is in agreement with the one used for estimating 
the magnetic field (v gaj = 1000 km s _1 , v gas = 905+}2j km s _1 : 
van Weeren et al. 2010; Stroe et al. 2014a). 

Here we discuss possible systematics on our estimate of pre 
shock temperature. The Mach number we determined is in prin¬ 
ciple also affected by any temperature gradient that was present 
before the shock passage. Any shock heating will appear as an 
excess above the temperature gradient that was already present. 
We can provide an approximate estimate of this effect by as¬ 
suming the gradient before shock passage followed the profile 
from Burns et al. (2010). Substituting an average temperature 
of <kT>-8 keV into Eq. 8 of Burns et al. (2010), the expected 
temperature gradient between the center of pre and post shock 
is 2.0 keV. So the temperature at the current post-shock loca¬ 
tion, but before shock passage was about 2.7 + 2.0 = 4.7 keV. 
After passage of the shock this increased to the currently mea¬ 
sured post-shock value of 8.5 keV. With a value of Ti=4.7 keV 
and f 2 —8.5 keV this corresponds to a Mach number of 1.8. 

Additional uncertainties arise from the fact that we assumed 
the temperature measurements are valid for the center of the 
bins in the above calculation. Our actual temperature estimate 


is emission weighted, resulting in a temperature bias in the high 
density region. Therefore, the position of the weighted temper¬ 
ature of the pre-shock bin could actually be located closer to 
the relic position (see fig. 6 in Hoshino et al. (2010)) and the 
shock heating is expected to be located at the edge of the relic. 
Due to this fact, the correction for the pre-shock temperature is 
likely somewhat smaller than what we assumed above. As an 
example, the expected temperature gradient between the edge 
of the relic and the center of the pre-shock region is ~0.8 keV, 
assuming the same Burns et al. (2010) profile. The expected pre¬ 
shock temperature precisely at the edge of the shock would then 
be 2.7+0.8=3.5 keV. With our measured post-shock temperature 
this would then corresponds to a Mach number of 2.3. We stress 
that the above calculations are just estimates as we lack informa¬ 
tion on the original temperature profile before the merger. The 
discussion in the following sections and conclusions are not af¬ 
fected by this systematics. 

3.3. Merger scenario of CIZA J2242.8+5301 

The temperature distribution of the ICM shown in Figure 3, 5, 
7 and the presence of the shock fronts suggest that a subcluster 
is colliding into the main body from the North-South direction. 
The peak of X-ray emission is close to the southern part of the 
cluster (Ogrean et al. 2013,2014), which might be infalling from 
the north direction. In this point of view, the northern shock cor¬ 
responds to the backward shock induced by the merging activ¬ 
ity. We do not see a significant excess in the temperature profile 
in the perpendicular direction to the merger axis, which is pre¬ 
dicted by several simulations (Akahori & Yoshikawa 2010; Mol- 
nar et al. 2013). Overall, our observations indicate that CIZA 
J2242.8+5301 is experiencing an almost head-on merger with 
near unity mass ratio and small impact parameter. 

Here we estimate the dynamical time scale of CIZA 
J2242.8+5301 by using its shock properties and radio infor¬ 
mation. Assuming that the shock propagation speed does 
not change from v s h ock ~ 2000 km/s and that the dis¬ 
tance from the north to the south relic is 2.5 Mpc, we esti¬ 
mate that CIZA J2242.8+5301 is observed approximately 2.5 
Mpc/(2xv J / !OCJ t) ~0.6 Gyr after core passage. Note that the ac¬ 
tual times scale is expected to be longer as the shock speed just 
after core passage is lower. Therefore in the end, the total travel 
time would be longer than 0.6 Gyr. Our results support pre¬ 
vious estimates based on hydrodynamic al simulation and radio 
observations (van Weeren et al. 2011; Stroe et al. 2014a). They 
suggested that CIZA J2242.8+5301 is a binary cluster merger 
with a mass ratio of 1:2, less than 10 degrees from the plane of 
the sky and impact parameter b < 400 kpc. The combination 
of radio and X-ray observations is powerful to constrain cluster 
merger events. For a more detailed understanding of the dynam¬ 
ical evolution of the merger, the mass distribution from a weak 
lensing observation and hydrodynamical simulations are useful 
(e.g., Dawson 2013). 

3.4. Mach number from X-ray and Radio observations 

Recent X-ray follow-up observations of the radio relics re¬ 
vealed a relation between the shock front and radio emis¬ 
sion (Finoguenov et al. 2010; Macario et al. 2011; Mazzotta 
et al. 2011; Ogrean et al. 2013; Ogrean & Briiggen 2013; Aka¬ 
matsu & Kawahara 2013). These results are consistent with the 
hydrodynamical simulations of merging clusters, which predict 
that merging events generate shock waves toward outskirts (Tak- 
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Fig. 10. The Mach number derived from radio (Mradio) plotted 
against that from the ICM temperature. The error bars are only sta¬ 
tistical and at the one sigma level. Red and blue crosses are this work, 
but with previous and new radio result (van Weeren et al. 2010; Stroe 
et al. 2014b). Grey crosses show the previous our results (Akamatsu & 
Kawahara 2013). Note: Stroe et al. (2014a) does not give information 
about the southern relic because the low flux. Therefore the result of 
the new radio analysis is only available for the northern relic. However, 
similar to the northern relic case, the Mach number of the southern radio 
relic is expected to drop. 

izawa & Mineshige 1998; Ricker & Sarazin 2001; Kang et al. 
2005). These shocks accelerate electrons up to relativistic ener¬ 
gies via diffusive shock acceleration (DSA) mechanism, which 
generates radio emission via synchrotron radiation. However, 
as mentioned in the introduction, the acceleration efficiency of 
DSA with low-M (< 10) shocks is known to be too low to 
take account for the radio flux. Therefore the main acceleration 
mechanism at cluster merger shocks is still unclear and until now 
the subject of discussion. For understanding of shock accelera¬ 
tion, it is crucial to know the particle acceleration efficiency and 
injection rate. 

Because of the presence of the two clear and giant radio 
relics, C1ZA J2242.8+5301 is one of the best targets to study 
shock acceleration in clusters. Especially, the northern relic 
shows an extremely narrow shape, whose width is only 50 kpc 
(Figure 1). van Weeren et al. (2010); Stroe et al. (2013) reported 
the spectral index map of the radio relics in C1ZA J2242.8+5301. 
At the edge of the northern radio relic, the spectral index is es¬ 
timated from the spectral index map and color-color diagram as 
a„ = -0.6 + 0.05, which correspond to the Mach number of 
M - 4.6 '((q. They also confirm the presence of a spectral index 
gradient across the relic, which is consistent with the DSA the¬ 
ory combined with Synchrotron aging. On the other hand, for 
the southern radio relic, it was not possible to estimate the injec¬ 
tion spectral index directly from spectral index maps due to poor 
signal-to-noise. Therefore, the estimated injection spectral index 
from the measured integrated spectral index a,,,, = -1.29 + 0.05 
as ai„j - -0.79 + 0.05, which corresponds to a Mach number of 
M = 2.8 + 0.2 (Stroe et al. 2013). We note that the estimated 
injection spectral index remains uncertain due to the assump¬ 


tions that have to be made if it is derived from the integrated 
spectrum (Kang 2015a,b). The new results of the radio anal¬ 
ysis (Stroe et al. 2014a) is slightly different from the previous 
result (a„ ~ -0.6: van Weeren et al. 2010; Stroe et al. 2013). 
The main difference is that Stroe et al. (2014a) used a largest 
radio data set and took into account spatial convolution effects. 
Because of the low flux of the southern radio relic, Stroe et al. 
(2014a) could not estimate the spectral index. Therefore the re¬ 
sult from the new radio analysis is only available for the northern 
relic. 

First order Fermi acceleration (under the assumptions of sta¬ 
tionary and continuous injection) gives relativistic electrons with 
a power-law spectrum n(E)dE ~ E p dE, where the power-law 
index is p — (C + 2)/(C - 1). The radio spectral slope just 
after the shock acceleration is a - —(p - l)/2. Using the mea¬ 
sured shock compression, this relation implies p n = 2.64 + 0.30, 
p s - 3.96 *qH and a„ = -0.81 + 0.15, a s = respec¬ 

tively. For the northern radio relic, the above value is consistent 
with the spectral index from the radio investigation (van Weeren 

et al. 2010; Stroe et al. 2014a, a -0.6 + 0.05, a -0.77). For 

the DSA theory, the measured injection index would require a 
shock compression of C„ ~ 2.84, which is again consistent with 
what is suggested by the temperature drop across the radio relic 
(Table. 5: C n ~ 2.8). 

Although it is likely to change similar to the northern ra¬ 
dio relic case, there is a difference between the X-ray and radio 
measurement of the southern radio relic. Using the measured in¬ 
tegrated spectral index and an assumption for the cooling effect, 
a similar trend has been seen in several radio relics and edges 
of the halo (Markevitch et al. 2005; Ogrean & Briiggen 2013; 
Ogrean et al. 2013; Akamatsu et al. 2013). There are several 
mechanism proposed to explain the particle acceleration in clus¬ 
ters (Markevitch et al. 2005; Kang et al. 2012; Pinzke et al. 2013; 
Skillman et al. 2013). A plausible possibility is re-acceleration 
of pre-existing non-thermal (low energy cosmic-ray) particles in 
the ICM. In this case, the pre-shock ICM already contains sev¬ 
eral type of non-thermal particles, which are generated by sev¬ 
eral mechanisms such as large scale accretion shock (Ryu et al. 
2003), turbulence acceleration (Brunetti 2011) and past activity 
of radio or starburst galaxies (see Briiggen et al. 2012, for a re¬ 
view). These low energy cosmic-ray particles have a long life 
time comparable to the ages of clusters or longer (Sarazin 1999; 
Pinzke et al. 2013). Recently, Kang et al. (2012) performed 
time-dependent simulations for diffusive shock acceleration for 
the case of the northern radio relic in CIZA J2242.8+5301. 
Compared with the observed synchrotron flux and spectral dis¬ 
tributions, they concluded that if a pre-existing population of 
low-energy cosmic-ray electrons exists, the radio relics can be 
explain by weak shocks with a Mach number M ~ 2-3. A 
similar conclusion is also reported by Pinzke et al. (2013). They 
found that contributions of fossil electrons can not be ignored for 
low Mach number shocks and their self-similar analytic model 
described the observed radio properties well. Another possibility 
is that of a varying Mach number along the shock. The ICM in 
the cluster outskirts is thought to be clumpy (Nagai & Lau 2011; 
Simionescu et al. 2011), the inhomogeneities in the ICM will 
thus result in small Mach number variations which could lead to 
a discrepancy between the radio and X-ray derived Mach num¬ 
bers as the shock acceleration efficiency (and radio relic bright¬ 
ness) scales non-linearly with Mach number (Hoeft & Briiggen 
2007). 

For CIZA J2242.8+5301, we measured the ICM temper¬ 
ature jump across the northern relic, which corresponds to a 
Mach number of Wf ~ 3. Indeed, the latest radio investiga- 
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tions by Stroe et al. (2014a) suggest that the previous radio 
spectral index might be biased and that the estimated spectral 
index will be a ~ -0.77, which correspond to a Mach num¬ 
ber M ~ 3.0. The new Mach number from radio observation 
agrees well with X-ray observations. In addition, the measured 
Mach number matches the prediction of the simulation by Kang 
et al. (2012). This supports the scenario of re-acceleration of the 
pre-existing non-thermal particles in the ICM. Upcoming LO- 
FAR low-frequency radio observations will shed new light on 
this problem. Such a high quality radio spectrum allows us to 
observe the injection spectral index directly. 

3.5. Magnetic field pressure at the northern radio relic 

As suggested by previous radio observations (van Weeren et al. 
2010), the northern relic has a relatively high magnetic field 
(B ~ 5 pG). Similar results are also reported for other radio 
relics (Nakazawa et al. 2009). Such a magnetic field strength 
is comparable to the expected value at the cluster central region 
from simulations of Faraday Rotation Measure in the Coma clus¬ 
ter (Bonafede et al. 2010). 

The thermal pressure of the ICM and the magnetic pres¬ 
sure of relativistic electrons is calculated as / J ,h - kT ■ n e and 
Pmagnetic - B 2 /( 87t), respectively. Assuming the ICM tempera¬ 
ture and electron density kT ~ 8.0 keV and n e ~ 5 x 10~ 4 cm~ 3 , 
we estimate the thermal pressure at the northern radio relic 
P t h ~6.7 eV/cm 3 , assuming that the mean molecular weight is 
0.6. Here we assume that all X-ray emission comes from the 
thermal component. With a magnetic field strength of B ~ 5 
pG. the energy density is estimated as P magn etic ~ 0.62 eV/cm 3 . 
Comparing both pressures, we find that the magnetic pressure 
reaches ~ 9 % of the thermal pressure. In the above estimate, 
we did not include the contribution of the energy density of the 
non-thermal electrons and turbulence. Therefore PsIPth ~ 0.09 
is a robust lower limit of the non-thermal pressure. 

3.6. Possibility of non-equilibrium effects 

Finally we discuss the possibility of non-equilibrium effects. 
Similar to supernova remnants, cluster merger shocks are ex¬ 
pected to be collisionless. In the shock in young SNR, there 
is a well known relation between the ratio of the ion and elec¬ 
tron temperature T e /Tt and shock velocity as T e ITj ~ v~ 2 lock (van 
Adelsberg et al. 2008). In the cluster case, non-equilibrium 
effects such as non-Maxwellian electron distributions (Kaastra 
et al. 2009) and non-equipartition of electrons and ions tem¬ 
perature (Takizawa 2005; Rudd & Nagai 2009; Akahori & 
Yoshikawa 2010) and non-equilibrium ionization are also ex¬ 
pected. 

Because of the different mass between the ions and the elec¬ 
trons, in general, shock heating is more effective for the ions. 
The electrons set the same temperature via Coulomb collisions. 
After that, electrons stabilize down to a Maxwellian distribution. 
Finally, via ion to ion collisions, the ionization balance in the 
ICM reaches equilibrium. The collisional time scale is given by 
n e ti e = 3 x 10 11 cm~ 3 s. Here, n e t le is the ionization parameter, 
which is an indicator of the ionization degree (Masai 1984). The 
region just behind the shock may have a higher ion temperature 
T, than the electron one T e , because it has not had enough time 
to reach equilibrium. We also tried to evaluate the ionization 
parameter n e t at the post-shock region using the NIE model in 
Xspec. Unfortunately Suzaku can not constraint well the n e t pa¬ 
rameter because of the low abundance and low statistics of Fe- 


K lines. Further studies should be possible by high-resolution 
spectroscopy of the ICM with the Japanese ASTRO-H (Taka- 
hashi et al. 2012) and the European ATHENA 10 (Nandra et al. 
2013) satellites. 

4. Summary 

CIZA J2242.8+5301 (z = 0.192) hosts two well defined giant 
radio relics. We find that the temperature profiles show a re¬ 
markable drop across both relics. Below is a summary of our 
main results: 

- The ICM temperature of the merging direction shows a flat 
profile with kT ~ 8 keV and shows a clear drops to 3-5 keV 
across the radio relics. The ICM temperature of the perpen¬ 
dicular direction to the merging axis shows different profiles. 
The temperature profile in the western direction agrees with 
the simulated temperature profile of Burns et al. (2010). 

- The significant drop in the ICM temperature indicates the 
presence of shock fronts. Based on the temperature drop, the 
estimated Mach number of the shocks are M„ = 2.7 l0 ' 3 , and 
M s = 1.7+ 0 - 4 , respectively. The shock velocities are also es¬ 
timated as v shockm ~ 2300 km/s and v s h oc k:s ~ 2000 km/s, re¬ 
spectively. This suggests that the merger happened ~0.6 Gyr 
ago, in agreement with estimates based on hydrodynamic al 
simulation and dynamical analysis (van Weeren et al. (2011); 
Stroe et al. (2014a) and Dawson et al. in prep.). 

- The observed radio spectral index of the northern relic (Stroe 
et al. 2014a) is consistent with what would be expected from 
diffusive shock acceleration in a shock with the observed 
Mach number and compression. 

- By combining our X-ray data with radio results, we es¬ 
timated the fraction of the magnetic pressure to the ther¬ 
mal pressure at the northern radio relic to be about ~ 0.62 
eV/cm 3 , which corresponds to 9% of the thermal pressure. 
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